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ABSTRACT 

We present a generalization of the multiphase chemical evolution model applied to a 
wide set of theoretical galaxies with different masses and evolutionary rates. This general- 
ized set of models has been computed using the so-called Universal Rotation Curve from 
Ifersic, Salucci & Steel ( 1996) to calculate the radial mass distribution of 44 theoreticaJ proto- 
galaxies. This distribution is a fundamental input which, besides its own effect on the galaxy 
evolution, defines the characteristic collapse time-scale or gas infall rate onto the disc. We 
have adopted 10 sets of values, between and 1, for the molecular cloud and star formation 
efficiencies, as corresponding to their probability nature, for each one of the radial distribu- 
tions of total mass. Thus, we have constructed a bi-parametric grid of models, depending on 
those efficiency sets and on the rotation velocity, whose results are valid in principle for any 
spiral or irregular galaxy. The model results provide the time evolution of different regions of 
the disc and the halo along galactocentric distance, measured by the gas (atomic and molec- 
ular) and stellar masses, the star formation rate and chemical abundances of 14 elements, for 
a total of 440 models. This grid may be used to estimate the evolution of a given galaxy for 
which only present time information - such as radial distributions of elemental abundances, 
gas densities and/or star formation, which are the usual observational constraints of chemical 
evolution models - is available. 

Key words: galaxies: abundances - galaxies: evolution- galaxies: spirals -galaxies: stellar 
content 



1 INTRODUCTION 

Chemical evol i jtion i nodels (CEM) jLvnden-Belllll975l: Finslevl 
Il980': Clavton'. ligSlfTgiSlSommer-Larsen & Yosh ii. 1989) were 
early developed to try to understand the origin of the radial gradi- 
ents of abundances, observed in our Galaxy (MWG). Most numer- 
ical models in the literature, including the multiphase model used 
in this work, explain the existence of this radial gradient by the 
combined effects of a star formation rate (SFR) and an infall of gas 
which vary with galactocentric radius in the Galaxy. 

A radial decrease of abundances has also been observed 
in most spiral galaxies I Henrv & Worthevl Il999h although the 
shape of the radial distribution changes from galaxy to galaxy. 
Among other global trends it is found that for isolated non- 
barred spirals the steepness of the radial gradient depends on 
morphological type, with later types showing steeper gradi- 
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ents jPiazl Il989h . with other general galaxy properties as sur- 
face brightness and neutral and molecu l ar gas fractions also 
playin g a role jVila-Costas & Edmundsl Il992t IZaritskv et all 
11994) . The radial gradient tends to be wiped out however for 
strongly barred galaxies which show flat abundance distribu- 
tions I Martin & Roy, 1994; Roy, 1996; Friedli, Benz, & Kennicutl 

r 

1994). Irregula rs galaxies also show uniform abundances t h rough- 
out (Rov. Bellev. Dutil. & Ma rtirl 11994 IWalsh & Ro^l 1 19971 
Kobulnickv, 1998; Mo lla & Ro4ll999h . 

The abundance gr adient pattern seems to show an on-ofF mode 
iEdmunds & RovllT993i) . being very steep for the latest spiral types 
and very flat for irregulars. All these considerations become clear 
when the gradient is measured in dex/kpc, but there are indications 
that suggest a gradient independent of galaxy type when it is mea- 
sured in dex/scale length (Diaz, 1989; Garnett, 1998). In order to 
analyze the behaviour of the radial distribution of abundances and 
the value of the radial gradient from a theoretical point of view a 
large number of models is necessary. Historically, CEM aiming to 
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reproduce radial abundance gradients iiave been, iiowever, applied 
only to the MWG. 

Actually, there is a lack of tools to determine the chemi- 
cal evolutionary state of a particular galaxy, besides our works 
applying the multiphase models to spiral galaxies. The recen t 
works bv lBoissier & Prantzo^ <200(i) : |Prantzos & BoissieJ <2000h 
are valid for galaxies other than the MWG. Their calculations 
use the angular momentum and rotation curves as model inputs 
keeping the star formation efficiency constant for all galaxies 
iBoissier. Boselli. Prantzos. et all 1200 ih . This technique may not 
be flexible enough to validate the models against observational 
data. In fact, a comparison to see if these models reproduce the 
observed abundance radial distributions of particular galaxies has 
not been done. It is always possible to extract some informa- 
tion by using evolutionary synthesis models in comparison with 
spectro-photometric observations. This method, very useful for the 
study of elliptical galaxies, does not result equally successful in the 
case of spiral galaxies due to the difficulty of measuring the spec- 
tral indices, exce pt for the bulges iOoudfrooii. Goreas. & Jablonkal 
1 19991 |Proctor& Sansom, 2003), from which ages and metallic- 
ities are obtained. Furthermore, even when these measurements 
are do ne with confidence iBeauchamD & Hardv, 1997; Molla et aLj 
1 19991) . in order to apply this technique to spiral galaxies, a com- 
bination of chemical evolution and evolutionary synthesis models 
is required to solve the uniqueness problem associated to the first 
ones and the age-metallicity degeneracy associated to the second 
ones iMplla & Hardy, 200J. 

At present, the available options are either to use the clas- 
sical closed box model or a Galactic Chemical Evolution (GCE) 
model. However, the closed box scenario is recognised to be 
inadequate to describe the evolution of most galaxies and in 
fact its application in many cases can yield misleading results 
JValle. Ferrini. Galli & Shorel l2002h . In particular, the fact of as- 
suming that a system has a constant total mass with a monoton- 
ically decreasing star formation according to a Schmidt law, pre- 
vents the reproduction of the observational characteristics of most 
galaxies. On the other hand, the evolution of a galaxy with present 
time properties different from the Milky Way will not necessar- 
ily be equal to that predicted by a GCE model. Realistic chemical 
evolution models adequate to describe different types of spiral and 
irregular galaxies are therefore clearly needed. 

The multi phase model, whose charac teristic s have been de- 
scribed in iFer rini. Matteucci, Pardi et alj il992h . has been ap- 
plied and checked against obser vational constraints not only 
for the Milky Way Galaxy tFerrini. Molla. Pardi et all 1 1994 
iMoUa & Ferrinil Il995h . as it is commonly done, but also for 
a sample of spiral galaxies (discs and bulges) of different 
morp h ological types and total rn asses I Molla, Ferrini & Diaz, 
1996; Molla. Hardy & Beauchanicl Il999i: iMolla & Rovt 119991; 
Molla, Ferrini & Gozzj, 2000). The observed radial distributions 



of gas, oxygen abundances and star formation rate have been re- 
produced rather successfully and the observed correlations be- 
tween abundan ce g radients and galax y characteristics are also re- 
produced <MoTlTet'ari 119961: IMoIM&RovL Il999bt) . This galaxy 
sample, which includes the best studied objects, is however small 
(only 11) and encompasses a restricted range of morphologies and 
masses. The application of the model can however be extended to a 
larger sample if an adequate parameter space is defined thus provid- 
ing the required chemical evolution of different types of galaxies. 

The model uses as input parameters the collapse time scale to 
form the disc, which depends on the total mass of the galaxy, and 
the efficiencies to form molecular clouds and stars which we as- 



sume different from galaxy to galaxy. The radial distributions of to- 
tal mass constitute the fundamental input of the multiphase model. 
They are easily computed when the rotation curves are available 
(Molla & Marquez, in preparation). If this is not the case, some 
assumptions are necessary. In this work, we have used the Uni- 
versal Rotation Curve from lPersic et al ]jl99dhereafter PSS96) to 
calculate a large number of mass radial distributions representing 
theoretical protogalaxies or initial structures which will evolve to 
form the observed spiral discs or irregulars. The total mass of each 
simulated galaxy, besides having its own effect on the galaxy evo- 
lution, defines the characteristic collapse time-scale or gas infall 
rate onto the disc. Regarding molecular cloud and star formation 
efficiencies, which will take values between and 1, we have cho- 
sen 10 different sets of values for each radial distribution of total 
mass. We have computed how the chemical evolution proceeds for 
galaxies defined by the different parameter combinations. The final 
bi-parametric grid consists of 440 models simulating galaxies of 44 
different total masses. 

This work represents an extension of our previous work that 
can help to understand the general trends observed in spiral and ir- 
regular galaxies concerning neutral and molecular gas distributions, 
abundance radial distributions etc. But, most importantly, by using 
these models we can predict the time evolution of a galaxy when 
only present time data are known. 

In Section 2 we summarize the general characteristics of the 
multiphase chemical evolution model and the strategy of its appli- 
cation to the objects of the grid. The results are presented in Section 
3 including the time evolution of galaxies and some radial distribu- 
tions for the present time. In section 4 we give the calibration of the 
grid models with the MWG and a restricted sample of well studied 
spiral galaxies and we discuss the model results in a global way. 
Finally, the conclusions of this work are presented in Section 5. 

The results obtained in this work will be available in elec- 
tronic form at CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5), via http://cdsweb.u-strasbg.fr/Abstract.html, or at 
http://wwwae.ciemat.es/~mercedes. 



2 THE MULTIPHASE MODEL 

The model used in this work is a generalization of that developed 
for the Solar Neighborhood in Ferrini et al. ( 1992) and later applied 
to the whole MWG l ,Fe rrini et al., 1994) and other spiral galaxies 
iMolla et allll996lll999lh ^ The enriched material proceeds from the 
restitution by dying stars, considering their nucleosynthesis, their 
initial mass function (IMF) ~ and hence the delayed restitution - 
and their final fate, via a quiet evolution, or Type I and II super- 
nova explosions. Most recent works support the idea that the IMF 
is practically universal in space a nd constant in time <Wvsd.ll997t 



We have adopted the IMF from 


Ferrini, Palla & Penco 1 1990), very 


similar to a Scalo's law 


Scalo 


1986) and in good agreement with 



the most recent data from iKroupa l2001i) . as can be seen in Fig.0 
The original model has been modified in order to use metal- 
licity dependent yields. Nucleosynthesis yields for massive stars 
have been taken from Wooslev & Weaver 1 1995). For low mass 
and intermediate stars we h ave used the set of yields from 
iGavilan. Buell & Moll^ i2004[l For type I supernova explosion re- 
leases, the model W 7 from iNomoto et al, 1 19841) . as revised by 
llwamoto et alj i 19991) has been taken. 
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Figure 1. The IMF from iFerrini. Palla & Peiicol ll99d) . solid li ne, com- 
pared to a Salpeter law, dotted line, and that corresponding to Kroupa 
I2OOII) . short-dashed line. Solid symbols correspond to .Scalg 1.1986.) . 
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Radial distributions: a) total masses M,o,, b) masses AM{R) in- 
our cylinders, for different values of A following the labels in the 



2.1 The rotation curves and mass radial distributions 

In this model each galaxy is described as a two-zone system with a 
halo and a disk. It is assumed that the halo has a total mass which 
is initially in gas phase. The total mass, M, and its radial distribu- 
tion, M(R), are calculated from the corresponding rotation curve 
derived from the Universal Rotation Curve of PSS96. These au- 
thors use a homogeneous sample of about 1100 optical and radio 
rotation curves to estimate their profile and amplitude which are 
analysed statistically. From this study they obtain an expression for 
the rotation velocity, V(R), as a function of the rotation velocity at 
the optical radius (the radius encompassing 83% of the total inte- 
grated light), Vopi, the galaxy radius normalised to the optical one, 
X = R/Ropt, and a parameter A = L/L,, which represents the ratio 
of the galaxy luminosity to that of the MWG, L, = 10'" '' Lq: 

V(R) = V,ptF{A) kms-' (1) 

where 

FiA) - {(0.72 . 0.441og.)^-,^^^ . ,.ee--^^—-—^ 

These authors also found that Vopt in kms"', depends on /I as : 

_ 200^» '" 

°'" ~ [0.80 -I- 0.49/og^ + (0.75e-o"*)/(0.47 + 2.25/10*)] 1/2 

The same occurs with the optical radius, which depends on 
luminosity through the expression: R^pt = ISCL/L,)"^ kpc, which 
we will also use. 

In Table □ we show the characteristics obtained with these 
equations for 44 different values of A. In column (1) we give 
the number of the radial distribution, defined by the value of A, 
given in column (2). The optical radius, Ropt, and the virial radius, 
ifgai = 14.8/l"'' '''/?opt, are given in columns (3) and (4) respectively. 
We have defined a characteristic radius for each galaxy, which we 
will use as our reference radius, as R^ = i?opt/2. This radius, in 



kpc, is given in column (5). Column (6) gives the rotation veloc- 
ity, in kms"', reached at a radius Rm = Ropt/l-45 kpc (see PSS96 
for details). The total mass of the galaxy, calculated with the clas- 
sical expression Mg^i = 2.3210'^/?g,iV'i„ (Leaueuxl fTiil . m units 
of 10' Mq, is given in column (7), and, finally, the characteristic 
collapse time scale, in Gyr, which will be described below, is given 
in column (8). 

Each galaxy is divided into concentric cylindrical regions 1 
kpc wide. From the corresponding rotation curve we calculate 
the radial distributions of total mass, M(R), with the expression 
M(R) = 232lO^RV(Rf. ' 

From these distributions of the total mass, we easily obtain the 
one corresponding to each of the cylinders, AM(R). Both distribu- 
tions M{R) and AM(R) are shown in Fig.|5| The total mass includes 
the dark matter component (DM) which, in principle, does not 
take part in the chemical evolution. However,the DM contribution 
seems to be negligible for the large massive galaxies, mostly in the 
regions where the chemical evolution is calculated. According to 
/2 "Palunas & Williams ( 2000); Sellwoo d & Kosows kv ( 2000,and ref- 
(2^rences therein) 75% of the spiral galaxies are well fitted without 
a dark matter halo and the failure to reproduce an other 20% is di- 
rectly related to the existence of non-axisymmetric structures (bars 
or strong spiral arms). 

2.2 The collapse time-scales 

The gas computed following the previous section collapses to fall 
onto the equatorial plane forming the disc as a secondary structure. 

' This equation is vaUd only for spherical distributions. The correspond- 
ing expression for cylindrical distributions differs from the spheric al one 
only by a geometrical factor, between and 1 1 Burstein & Rubl^, Il985t 
ICampos-Aguilar. Prieto. & Garci'alll993ft . This difference is well inside the 
possible uncertainties of the actual rotation curves compared with the uni- 
versal rotation curve from PSS96. 
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Table 1. Galaxy Characteristics dependent on the Total Mass. 
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The gas infall from the halo is parametrised by fgn, where gu is 
the gas mass of the halo and / is the infall rate, that is, the inverse 
of the collapse time scale r. 

The collapse time-scale for each galaxy depends 
on its total mass through the expression: t oc M'^^'^Tg 
iGail agher. Hunter & Tutukov. 1984), where Mc, is the total 
mass of the galaxy in 10' Mq (column 7 of Table and Tg is 
its age, assumed to be 13.2 Gyr in all cases. We calculate a 
characteristic collapse time scale for each galaxy from the ratio 
of its total mass, A/g g£y/, and the MWG, A/q.mivg' 

Tc = TcMWciMg^gal / Mg^MWcT^'^ (4) 

where t^mwo is the collapse time scale at the characteristic 
radial region for the MWG, = Ropil'2- = 6.5 kpc. This col- 
lapse time scale (~ 4 Gyr) was determined from the corresponding 




100 200 300 

Vmax (km/s) 




5 10 15 20 25 



Radius (kpc) 

Figure 3. a) Characteristic collapse timescale r^. for each galaxy accord- 
ing to the maximum rotation velocity, b) Radial distribution of the collapse 
times scales TcoIi{R)- Lines have the same meaning than in Fig.|2|The hori- 
zontal solid line represents the assumed age of galaxies, 13.2 Gyr. 



value for the Solar Neighborhood, Tq ~ 8 Gyr. This value for Tq is 
very similar to that found in other standard galactic chemical evo- 
lution models jP ortinari. C hiosi. & Bressan [ ll998l:IChiaDDini et al.L 
ll99gtlChang' et al... .1999i) and it is constrained by a large number 
of data, such as the [O/Fe] vs [Fe/H] relation for stars i n the halo 
and the disc, the present infall rate, (0.7 Mapc^^Gyr"' , see'Mirabel 
1989«), or th e G-dwarf me tallicity distribution jperrini et al., 199^ 
IPardi&Fer 'rini. 1994; Kot oneva et"ail l2002h . all of them well re- 
produced in our mo del with this long scale to form the disc (see 
iGavilan et all Eo04l) . We would like to emphasize that the char- 
acteristic collapse time-scale computed with equation |4]is not a 
free-fall time ( < 1 Gyr for the MWG), being much longer than 
that since it has been calculated through the calibration with the 
Solar Neighbourhood collapse time scale Tq. 

The mass which does not fall onto the disc will remain in 
the halo, and yields a ratio M 1,^10 /M^n,^ for the baryonic compo- 
nent which is also in agreement with observations. The relative 
nor malization of the halo, thick and thin disc surface mass densi- 
ties teandagelfigSTl) gives an approximated proportion of 1 : 22: 200, 
which implies that the halo surface mass density must be about 
1/100 that of the disc component, also in agreement with the ratio 
obtained by the multiphase model. 

Fig. a) shows the characteristic time-scales vs the rotation 
velocity Vop, for our models. Solid dots represent the values used 
in our previous m odels for individual spiral galaxies 1 Ferrini et akl 
ll994HlMolla etarill99lll999ll . 

An important consequence of the hypothesis linking the col- 
lapse time scale with the total mass, is that low mass galax- 
ies take more time to form their discs, in apparent contra- 
diction with the standard hierarchical picture of galaxy for- 
mation. This characteristic is, however, essential to repro- 
duce most observational features of spiral and irregular galax- 
ies (see also Boissier & Prantzos, 2000). In fact, recent self- 
consistent hydrodynamical simulations in the context of a 
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cosmological model ^omm guez-Tenreiro. & Sernal l2004t 

iDomfnguez-Tenreiro. Saiz' & SemalllOoir show that a large pro- 
portion of massive objects are formed at early times (high redshift) 
while the formation of less massive ones is more extended in time, 
thus s imulating a modem version of the mon olithic collapse sce- 
nario iEggen. Lvnde n-Bell. & Sandagel ll962l ELS). Furthermore, 
this assumpti on is supported by various sets of observations 
fcimatti et al.', 2002; Jimenez et al., 2004; "Heavens et al.', '2004 
[ Olaze brook et al.. 2004; Cimatti et al., 2004; McCarthv et all 
l2004h which seem to demonstrate that a large proportion (~ 80 
%) of the massive galaxies formed their stars at z > 1 while a much 
smaller proportion of the less massive ones have converted their 
baryon mass into stars at that redshift. 

It is evident that the collapse time-scale must vary with galac- 
tocentric radius. Assuming that the total mass surface density fol- 
lows the surface brightness exponential shape, the required collapse 
time-scale should also depend exponentially on radius, with a scale 
length Aq cc Rq (where Rp is the corresponding scale length of the 
surface brightness radial distribution). Thus, we assume: 



Tcoiii^) = Tc exp ((R - R,)/iD) (5) 

In principle, w e might expect that X[, = 2Rd- However, in 
lFerrinietal] jl994^ we have shown that such large scale length for 
the infall rate produces a final radial variation for the elemental 
abundances in disa greement with the observ ed distributions in spi- 
ral galaxies, (see iPortinari & Chiosil [l99^ for a wide discussion 
about the effect of this parameter on the radial distribution of abun- 
dances). On the other hand, we must bear in mind that the surface 
brightness distribution is the final result of the combination of both 
the collapse and the star formation processes, and therefore the col- 
lapse time-scale may have in principle a different dependence on 
radius than the surface brightness itself. For the sake of simplicity, 
we assume a scale length A[, = OASR^,,,, corresponding to half the 
scale length of the exponential disc Rd = Roptl^-'2- given by PSS96. 
The value of Ad de creases with the mass of the galaxy, in agreement 
with observations I'Simien & de Vau couleurslll986l ; lGuzman et 2A\. 
Ll996; Graham & de Block, 2001,see their Fig. 4). 

The radial dependence of the infall rate is not imposed a priori 
in our scenario, but it is consequence of the gravitational law and 
the total mass distribution in the protogalaxy. The physical mean- 
ing is clear; galaxies begin to form their inner regions before the 
outer ones in a classical inside-out scheme. This halo-disc connec- 
tion is crucial for the understanding of the evolution of a galaxy 
from early times, the inside-out scenario being essential to repro- 
duce the radial grad i ent o f abundances IPortinari & Chiosi, 1993; 
Poissier & Prantzos, 2000). In fact, in a chemo-dynamical model 
cSamland. Hensler & Theis..l997.^ . this scenario is produced natu- 
rally. 

The radial variation of the collapse time-scale for each galaxy, 
calculated with equation |5| is shown in Fig.|3ji), where we also 
draw a solid line at 13 Gyr, the assumed age of galaxies. If the 
collapse time-scale is larger than this value, there is not enough 
time for all the gas to fall onto the disc; only a small part of it has 
moved from the halo to the equatorial plane and the disc formation 
is not yet complete. This could correspond to the situation observed 
bv lSancisi. Fraternali. Oosterloo. & van M oorse (2001), who have 
found an extended component of Hi, different from the cold disc, 
located in the halo, rotating more slowly than the disc and with 
radial inward motion. 



2.3 The star and cloud formation efficiencies 

The model computes the time evolution of each population which 
inhabits the galaxy. In the various regions of the disc or bulge and in 
the halo, which are treated separately, we allow for different phases 
of matter aggregation; diffuse gas (g), clouds (c, except in the halo), 
low-mass stars (si,m < 4Mq), high-mass stars {S2,m ^ 4Mo), and 
stellar remnants ^. The value for stellar mass range division is re- 
lated to nucleosynthesis prescriptions; stars with masses lower than 
4 Mq only produce light elements, and do not contribute to the in- 
terstellar medium enrichment.^ 

The mass in the different phases of each region changes by the 
following conversion processes, related to the star formation and 
death; 

(i) Star formation by the gas-spontaneous fragmentation in the 
halo 

(ii) Cloud formation in the disc from diffuse gas 

(iii) Star formation in the disc from cloud-cloud collisions 

(iv) Induced Star formation in the disc via massive star-cloud 
interactions 

(v) Diffuse gas restitution from these cloud and star formation 
processes 

In the halo, devoid of molecular clouds, the star formation fol- 
lows a Schmidt law for the diffuse gas with a power n = 1.5 
and a proportionality factor K. In the disk, stars form in two steps; 
first, molecular clouds, Co, form out of the diffuse gas, go, also 
by a Schmidt law with n = 1.5 and a proportionality factor called 
yu. Then, cloud-cloud collisions produce stars by a spontaneous pro- 
cess at a rate proportional to a parameter H. Moreover, a stimulated 
star formation process, proportional to a parameter a, is assumed. 

An advantage of using our model is that it includes a more re- 
alistic star formation than classic chemical galactic evolution ones 
in which the SFR prescriptions are based on a Schmidt law, depend- 
ing on the total gas surface density. Instead, the multiphase model 
assumes a star formation which takes place in two-steps; first, the 
formation of molecular clouds; then the formation of stars. This 
simulates a power law for the gas density with an exponent n > I 
and with a threshold gas density as shown by iKennicuttI ^98^, 
and, more importantly, it allows the calculation of the two differ- 
ent gas phases present in the interstellar medium. In fact, the actual 
process of star formation, born by observations ( Klessen, 2001), is 
closer to our scenario with stars forming in regions where there are 
molecular clouds, than to the classical Schmidt law which depends 
only on the total gas density. 

Our assumed SFR implies that some feedback mechanisms 
are included naturally and are sufficient to simulate the actually ob- 
served process of creation of stars from the interstellar medium. 
The formed massive stars induce the creation of new ones. But, 
at the same time, these star formation processes also may destroy 
the diffuse or molecular clouds, thus preventing the total conver- 
sion of gas into stars and ejecting more gas once again into the 
ISM -point (v)-. In particular, massive stars destroy the molecular 
clouds that surround them, due to the sensitivity of molecular cloud 
condensation to the UV radiation (Parravano, 1990). This mecha- 
nism restores gas to the ISM, thus decreasing the star formation. 

^ All these quantities represent the total mass of each phase, included in 
the code in units of IO'Mq 

^ Also, this splitting of stars into two groups, less and more massive than 4 
Mq, allows a very easy comparison of our resulting metallicity distribution 
with the observed one, based on G-dwarf low mass stars. 
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Both regulating process are included in our model, although nei- 
ther heating or cooling mechanisms for the cloud components are 
included explicitly in our code. 

The complete set of equations is given in lMolla et alj il996h 
for each radial region. We summarize here only those related to the 
star formation processes in the halo, f^, and in the disk, 



dcp 
dt 
I'd 



flg'lf - aCDSD.2 - Hc% 

Hci, + acosoa 



(6) 
(7) 
(8) 



where H and D indicate halo and disk. K, ji, H and a , besides the 
previously described t are the parameters of the model. 

In the classical method of application of a chemical evolution 
model to a given region or galaxy, the input parameters are con- 
sidered as free and chosen as the best ones in order to reproduce 
the selected observational constraints of the galaxy. In our models, 
however, not all these input parameters can be considered as free. 
The parameter / is defined by the total mass radial distribution, 
as already explained. Regarding K, ja, H and a, we have tried to 
reduce to a minimum their degree of freedom. Their radial depen- 
dence may be estimated as iFerrini et alill994l see): 



K = 6^(G/\/«)'/2 
H = en(ctelVD) 



(9) 
(10) 
(11) 
(12) 



where G is the universal gravitational constant, Vh and Vd are the 
halo and disk volume of each radial region, is the average cloud 
density and < m.„ > is the average mass of massive stars. The pro- 
portionality factors are called efficiencies that represent probabili- 
ties associated with these processes. 

In this way, the free parameters K, ji, H and a, proportionality 
constants in our star and cloud formation laws, but variable for each 
radial region, are computed through the efficiencies, which repre- 
sent the efficiencies of star formation in the halo, e^, cloud for- 
mation, e^, cloud — cloud collision, eu, and interaction of massive 
stars with clouds, e„, in the disc. Only a given efficiency for each 
process must be selected for the whole galaxy, although the original 
parameters, and its corresponding processes, maintain their radial 
dependence. 

The term associated to the induced star formation describes a 
local process and, as a result, its coefficient e„ is considered inde- 
pendent of the galaxy modelled and the location in it. The term 
is also assumed constant for all haloes. Therefore, for all our 440 
models, both efficiencies take the same values already used in our 
previous model for the MWG. Only the other two efficiencies 
and Eh are allowed to vary for each galaxy, being characteristic for 
each one of them. Within each galaxy eh and are independent of 
the position. 

The fact that galaxies with the same gravitational potential or 
mass but different morphological type or appearance exist, implies 
that the evolution of a galaxy does not depend solely on gravitation, 
even if this may be the most important factor, but also on certain dy- 
namical conditions. These conditions cannot be taken into account. 

The volume of the disc is calculated with a scale height of 0.2 kpc for 
all galaxies; the volume of the halo for each concentric region is computed 
through the expression: Vi,;,io(R) = 2Rhaio ^(1 - T^f- 



obviously, in a simple chemical evolution model, but may change 
the evolution of a galaxy, (mostly the star formation rate through 
temperature variations). We may consider them as included in our 
efficiencies to form molecular clouds and stars, 6^ and eh, which 
are allowed to change from one galaxy to another. 

We assume that and vary between and 1, according 
with its efficiency meaning. In principle both parameters should 
be allowed to vary independently from each other, what would in- 
crease very much the number of models to be calculated. Further- 
more, some of those models will not be physically possible. On the 
other hand, on the basis of our previous works, a trend between 
both efficiencies seems to exist, increasing or decreasing together 
for a given galaxy respect to the values appropriate for the MWG 
model. We have therefore studied the observational data related to 
molecular cloud and star formation from the available gas in order 
to check if a correlation between these two processes exists. 

To this aim we have used the data from Youna et al. jl996h 
which correspond to a large sample of galaxies and refer to atomic 
and molecular gas masses and to luminosities in the IR band and in 
Ho- emission, two well known indicators of the star formation rate. 
With them, we have tried to establish if some correlation appears 
among our efficiencies, and e^- Since the modelled star forma- 
tion rate have two steps to form stars, we will analyse the transfor- 
mation of diffuse gas in molecular gas and of the conversion of the 
molecular clouds into stars. 

From equations (10) and (11), our parameters and e^, ef- 
ficiencies of transforming the diffuse gas in molecular clouds and 
forming stars out of them, can be written as: 



067 
, V 



0.3 



while from equations (7) and (8) we have: 
(dcldt) + SFR 



H 



SFR 



(13) 
(14) 

(15) 
(16) 



In these expressions the subscript D has been omitted, since in 
all cases we refer to the disc 

We may approximate ^ ~ j;, where At is the mean time nec- 
essary to transform the diffuse gas in a molecular cloud, in units of 
lO^yr: 



{c/At) + SFR 



SFR V 
~03 



0.67 



(17) 
(18) 



where g and c are the total mass in diffuse and molecular gas in 
units of 10' Mq, and SFR is the star formation rate in units of 
10''Mo/10''3'r. Therefore: 

log6^ = log((c/A0 + 5F^)- 1.51ogg + 0.51ogV + 0.18(19) 
log€n = log 5 - 2 log c + log V -I- 0.5 (20) 

SFR is u sually estimated from the Ha luminosity 
jYoung et all[l99Q^ : 

' Actually, T = /ic~ + acs2 but the second term is much smaller than the 
first and therefore we can approximate H ~ \ 
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Figure 4. The ratio Me^ //"£// as a function of tlie morpliological type T. 
Solid dots are the averaged value for each morphological type. The soHd 
line is the mean value estimated from these bins. 



S FRiW'' Mq/ 10' yr) = L„„3.410-"' (21) 

and the volume of the disc is 
V = O.lnR^ (22) 

where R is in kpc and Lh„ is in L©. We have all 
data to estimate these efficiencies, except At. Recent estimates 
(Prinde. Allen. & Lubow, 2001; Bergin et al., 2004) for this cloud 
accumulation time scale give values several times the gravitational 
contraction scale (which is ~ 1-4 Myr), that is At = 10 - 20 
Myr, and probably smaller than 50 Myr. The probable range is 

[3.10* 8.10'] yr. We have given three possible values: 1.10*, 

5.10' and 5.10* yr in order to take into account other possible 
slower modes (Palla, 2004). 

In Fig.0we represent the ratio Ine^/lneu computed for At = 
5. 10' vr as a function of galaxy morphological type. Single galaxies 
are represented as + while solid dots correspond to average values 
obtained by binning them into 11 types. A constant ratio for all 
types is consistent with the data points, shown by the solid line and 
giving a mean value < >~ 0.34(±0.02) with;^'^ = 3.6. * Taking 
into account the large range of variation of and (more than 5 
orders of magnitude) for the whole set of data, it is surprising that 
the Ine^/In^u ratio takes values around 0.35 for all galaxies. For 
other values of At, this value changes but remains constant. They 
are shown in the graph as dotted lines at values 0.00 (At = 1 Myr) 
and 0.80 (At = 50 Myr). Taking this into account,we have assumed 
a ratio ln£^/ln£u= 0.4 for our computed models. 

We have computed 10 models for each mass radial distribu- 
tion, allowing to take values between to 1 as given in Table|2| 
For each one of these efficiencies, 6/^ has been fixed according to 
the ratio < Z«e^//ne// >~ 0.4 and their values are also given in the 

* The increasing ratio for increasing type T is only apparent due to the 
reduced number of galaxies within the bins in types from T = 8 to 10. 



Table 2. Efficiencies chosen for our models. 



N 




e« 


1 


0.95 


0.88 


2 


0.80 


0.57 


3 


0.65 


0.34 


4 


0.45 


0.14 


5 


0.30 


0.05 


6 


0.15 


l.Oe-2 


7 


0.075 


1.5e-3 


8 


0.037 


2.6e-4 


9 


0.017 


3.7e-5 


10 


0.007 


4.0e-6 



table. Each set of efficiencies has been labelled by a number N from 
1 to 10. 

Summarizing, only the characteristic collapse time-scale, de- 
pending on the total mass, and the set of efficiencies, denoted by 
number N, are varied from model to model. We thus obtain 440 
different models which represent all possible combinations of the 
collapse time-scale with the values of N. 

3 RESULTS 

3.1 The time evolution 

The results corresponding to the mass of each region and phase, the 
star formation rate and the supernova rates, for the 440 computed 
models are shown in Tables |3| and |4| The elemental abundances 
for the discs are shown in Table |5| Here we only give, as an ex- 
ample, the results of the model corresponding to radial distribution 
number 22 and N = 5, with a rotation velocity of ~ 140 km.s"', for 
the first and last time steps of the evolution. 

The whole tables with the complete time evolu- 
tion from to 13 Gyr, with a time step of 0.5 Gyr, for 
the whole set of models will be available in electronic 
form at CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/Abstract.html, 
or http:/wwwae.ciemat.es/~mercedes/grid. 

In Table|21we list in column (1) the time, in Gyr and in column 

(2) , the galactocentric distance, in kpc, in Column (2). Columns (3) 
to (9) give the mass, in units of 10' Mq, in each region and phase: 
column (3) the total mass in each region; column (4) the mass of the 
disc region; column (5) the mass in the difi^use gas phase;column 
(6) the molecular gas; column (7) the mass in low and intermedi- 
ate mass stars; column (8) the mass in massive stars and, finally, 
column (9) the mass in remnants. 

Table|4|gives, for each time step in Gyr and radial distance in 
kpc , listed in columns (1), and (2) respectively , the star formation 
rate, in units of Moyr"', in the disc and the halo regions in columns 

(3) and (4); the supernova rates, for Types la, lb and and II, for the 
disc in columns (5), (6) and (7) and for the halo in columns (8),(9) 
and (10) all of them in units of 100 yr"' . 

The abundances in the disc for 14 elements are given in Ta- 
ble|5|for each time step in Gyr and radial distance in kpc (columns 
(1), and (2) ).The abundances by mass of: H, D, ^He, '♦He, '^C, '""C, 
N, O, Ne, Mg, Si, S, Ca, and Fe are given in columns (3) to (16). 

The time evolution of several models is shown in the next fig- 
ures. In each one, 4 panels are shown, corresponding to 4 different 
maximum rotation velocities and/or radial distribution of masses. 
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Table 3. Model Results corresponding to masses in each region and phase. 



Time R Mtot Mdisc Mgas(HI) Mgas(H2) Mstars(M< 4Mq) Mstars(M^ 4Mo) Mremnants 

(Gyr) (kpc) (lO'Mo) (IO'Mo) (ICMo) (ICMq) (IO'Mo) (IO'Mo) (ICMq) 



0.1 


14. 


0.69E+01 


0.47E-04 


0.47E-04 


0.89E-08 


0.40E-17 


0.31E- 


-18 


0.21E-19 


0.1 


12. 


0.69E+01 


0.22E-03 


0.22E-03 


O.llE-06 


0.65E-15 


0.49E- 


-16 


0.34E-17 


0.1 


10. 


0.68E+01 


O.lOE-02 


O.lOE-02 


0.14E-05 


0.13E-12 


0.93E- 


-14 


0.65E-15 


0.1 


8. 


0.67E+01 


0.48E-02 


0.47E-02 


0.20E-04 


0.28E-10 


0.21E- 


-11 


0.14E-12 


0.1 


6. 


0.66E+01 


0.22E-01 


0.22E-01 


0.31E-03 


0.75E-08 


0.56E- 


-09 


0.39E-10 


0.1 


4. 


0.58E+01 


0.91E-01 


0.87E-01 


0.47E-02 


0.21E-05 


0.16E- 


-06 


O.llE-07 


0.1 


2. 


0.29E+01 


0.21E+00 


0.17E+00 


0.39E-01 


0.23E-03 


0.17E- 


-04 


0.13E-05 



13.2 


14. 


0.69E4-01 


0.58E-02 


0.46E-02 


0.13E-02 


O.lOE-04 


0.22E- 


-07 


0.70E-06 


13.2 


12. 


0.69E+01 


0.27E-01 


0.16E-01 


O.lOE-01 


0.89E-03 


0.16E- 


-05 


0.67E-04 


13.2 


10. 


0.68E+01 


0.13E+00 


0.54E-01 


0.43E-01 


0.28E-01 


0.31E- 


-04 


0.25E-02 


13.2 


8. 


0.67E4-01 


0.57E+00 


0.14E+00 


O.IIE-fOO 


0.29E-F00 


0.21E- 


-03 


0.32E-01 


13.2 


6. 


0.66E4-01 


0.23E+01 


0.26E-I-00 


0.19E4-00 


0.16E+01 


0.79E- 


-03 


0.21E+00 


13.2 


4. 


0.58E-t-01 


0.49E+01 


0.18E+00 


0.18E+00 


0.39E+01 


0.88E- 


-03 


0.65E+00 


13.2 


2. 


0.29E-F01 


0.29E+01 


0.19E-01 


0.51E-01 


0.23E-I-01 


0.95E- 


-04 


0.50E+00 



We have selected the values corresponding to A =0.03, 0.19, 1.0 
and 2.5 corresponding to galaxies with rotation velocities of 48, 
100, 200 and 290 kms"', respectively, and representing typical ex- 
amples of spiral and/or irregular galaxies. For each panel we show 
the results for the 10 selected values of the efficiencies (or equiv- 
alently 10 rates of evolution) from N = 1, corresponding to the 
highest efficiency values and hence the most evolved models, to 
N = 10, with the smallest values and thus the least evolved ones. 

We have computed models for galactocentric radii up to 
l.SRopt with a step which depends on total galactic, being larger 
(up to 4 kpc) for the most massive modelled galaxies and smaller 
(only IKpc) for the lowest mass ones. Results are shown only for 
the region located at the galactocentric distance closest to the char- 
acteristic radius Rc defined in TableQ 

Fig. m shows the time evolution of the diffuse gas density. In 
all cases, the diffuse gas surface density is seen to increase rapidly 
at early times and then declines slowly. The first abrupt increase 
is a consequence of the infall rate of gas from the halo and hence 
does not depend on the model efficiencies. The increase is faster for 
more massive galaxies. The later decline is due to gas consumption 
in the process of molecular gas formation and hence depends on the 
efficiency e^, producing lower densities for the galaxies with higher 
efficiencies. 

Fig|6| shows the time evolution of the molecular gas phase 
whose formation is delayed respect to that of the diffuse gas. It 
is clear that the maximum density is reached later than that cor- 
responding to the diffuse gas density. For instance, for N = 5, 
& = 0.1 Gyr for Vrot = 290 and 200 km s"', while it is 0.3 Gyr for 
Vrot = 100 km s"' and reaches 6t = 1.2 Gyr when Vrot = 48 km 

Fig.Qshows the star formation rate history for the same cases 
as before. These histories result extraordinary different, even for 
equal efficiencies. Taking into account that the radial region shown 
is equivalent in all galaxies, it indicates that the primary agent driv- 
ing the time evolution of the SFR is the galaxy mass through the 
collapse time scale. On the other hand, the resulting star forma- 
tion history in most models, is smooth, showing values larger than 
IO^(MoGyr"'pc"^) only at the early evolutionary phases of the most 
massive galaxies with very high efficiencies. 

The time evolution of the oxygen abundance relative to hydro- 



gen, expressed as 12+log(0/H), for the models described above is 
shown in Fig|H| Present time oxygen abundances look very similar 
for all the galaxies with N ^ 5. The models with low efficiencies, 
A' > 7, show very low abundances although always larger than 
12+log(0/H)= 7, about the lowest observed oxygen abundance in 
galaxies. 

Fig|9| shows the evolution of [O/Fe] when [Fe/H] increases. 
The usual relation for the MWG can be taken to be represented by 
the N = 4 line from panel c), showing almost a plateau for low 
metallicities ([Fe/H] < 1.5) and then a decline toward the solar 
value. The plateau appears in the models c orresponding to massiv e 
massive galaxies with high efficiencies i Matteucci et al.L Il998h . 
The less massive galaxies show a slow decline without plateau for 
low values of N, and a steeper decrease for higher values. Thus, we 
may expect tha t [O/Fe] in the less evolved galaxies will soon reach 
the solar value jMatteucciri992) . 

3.2 Present day radial distributions 

We now analyse the present time results obtained with this bi- 
parametric grid of models. In order to do this, we start presenting 
only the present time radial distributions of gas, oxygen abundances 
and star formation rate which we will analyse in the following sub- 
sections. The radial distribution of gas and elemental abundances 
provide two basic model constraints since they are easily extracted 
from observations. In addition, in many cases, the star formation 
surface density and/or the surface brightness can also be derived. 

3.2.1 Radial distributions of diffuse gas 

The radial distribution of atomic gas surface density for the present 
time is shown in Fig. 1101 We can see that the atomic gas sur- 
face density show s a maximum somewhere along the disc, as it is 
usually observed (Cavatte et al.'.'l99(t lRhee & van AIbad3.ll996t 
[Sroeils & Rhee, 1997, and see Fig |17t . The value of this maximum 
depends on N: the models with the highest efficiencies have smaller 
gas quantities and their maximum densities are around 3-4 Mq/ pc^. 
For intermediate efficiencies (4 < W < 7), these maximum values 
rise to ~ 5 - 8 Mo/pc^. For all efficiencies the radial distributions 
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Figure 6. Time evolution of the logarithmic surface density of the molecular gas of the radial region defined by for 4 diff'erent mass distributions following 
the Vrot's values from each panel. In each one of them 10 models are represented, for the 10 different efficiencies, following labels in panel a). 
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Figure 7. Time evolution of the logaritiimic star formation rate surface density of the radial region defined by Rc for 4 different mass distributions following 
the Vrot's values from each panel. In each one of them 10 models are represented, for the 10 different efficiencies, following labels in panel d). 
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Figure 8. Time evolution of the oxygen abundance 12 + log{0/H) of the radial region defined by Rc for 4 different mass distributions following the Vrot's 
values from each panel. In each one of them 10 models are represented, for the 10 different efficiencies, following labels in panel d) of the previous figure. 
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Table 4. Model Results: Star Formation Histories and Supernova Rates. 



Time R SFR(disc) SFR(halo) SNJa (disc) SN-lb (disc) SN-11 (disc) SN-la (halo) SN-lb (halo) SN-II(halo) 

(Gyr) (kpc) (Moyr"') (Moyr^') (lOOyr-') (lOOyr-') (lOOyr-') (lOOyr"') (100}^-') (lOOyr-') 



0.1 


14. 


0.26E-15 


0.51E-01 


0.17E-17 


0.14E-18 


0.42E-14 


0.45E-02 


O.UE-Ol 


0.15E+01 


0.1 


12. 


0.43E-13 


0.52E-01 


0.27E-15 


0.14E-16 


0.69E-12 


0.46E-02 


O.UE-Ol 


0.16E+01 


0.1 


10. 


0.82E-11 


0.54E-01 


0.52E-13 


0.25E-14 


0.13E-09 


0.47E-02 


0.12E-01 


0.16E+01 


0.1 


8. 


0.18E-08 


0.57E-01 


0.12E-10 


0.55E-12 


0.30E-07 


0.51E-02 


0.12E-01 


0.17E+01 


0.1 


6. 


().49E-06 


0.63E-01 


0.31E-08 


0.15E-09 


0.80E-05 


0.55E-02 


0.14E-01 


0.19E+01 


0.1 


4. 


0.13E-03 


0.61E-01 


0.89E-06 


().44E-()7 


0.22E-02 


0.55E-02 


0.14E-01 


0.18E+01 


0.1 


2. 


0.14E-01 


0.28E-01 


O.llE-03 


().60E-05 


().24E+00 


0.26E-02 


().65E-()2 


0.84E+00 


13.2 


14. 


0.52E-05 


0.45E-01 


0.13E-05 


0.56E-05 


0.16E-03 


0.18E-01 


0.54E-01 


0.14E+01 


13.2 


12. 


0.38E-03 


0.46E-01 


O.lOE-03 


0.42E-03 


O.llE-01 


0.18E-01 


0.54E-01 


0.14E+01 


13.2 


10. 


0.74E-02 


0.46E-01 


0.23E-02 


0.86E-02 


0.22E+00 


0.18E-01 


0.55E-01 


0.14E+01 


13.2 


8. 


0.50E-01 


0.44E-01 


0.18E-01 


0.59E-01 


0.15E+01 


0.18E-01 


0.53E-01 


0.13E+01 


13.2 


6. 


0.19E+00 


0.29E-01 


0.73E-01 


0.22E+00 


0.56E+01 


0.13E-01 


0.35E-01 


0.86E+00 


13.2 


4. 


0.21E+00 


0.25E-02 


O.llE+00 


0.25E+00 


0.62E+01 


0.24E^2 


0.31E-02 


0.74E-01 


13.2 


2. 


0.22E-01 


0.17E-07 


0.29E-01 


0.28E-01 


0.68E+00 


0.14E-03 


0.27E-07 


0.51E-06 




Figure 9. The evolution of the relative abundance [0/Fe] vs [Fe/H] for the radial region defined by Rc for 4 different mass distributions following the Vrot's 
values from each panel. In each one of them 10 models are represented, for the 10 different efficiencies, following labels in panel a) 



are very similar independently of their galactic mass, except for 
those corresponding to /I < 0.03 (Vrot = 48kms~') which show 
much lower densities, except in the central region. 

In fact, a characteristic shown by all distributions is the sim- 
ilarity among models of the same efficiencies but different total 
mass, only scaled by their different optical radius. With the excep- 
tion of the A = 0.03 model, all the others show, for a same N, 
differences small enough as to simulate a dispersion of the data. 



For the less evolved theoretical galaxies (N > 7), the neu- 
tral gas distribution display a clear dependence on galactic mass. 
The maximum density values are always large, as expected, due 
to the small efficiencies to form molecular clouds, which do not 
allow the consumption of the diffuse gas, but this maximum den- 
sity increases from ~ ISMg/pc^ for A = 0.03 (Vrot = 48kms"') 
to ~ 30 - 40Ms/pc^ for A = 1.5 (Vrot = 248kms-') and up to 
~ lOOMo/pc^ for A = 2.5 (Vrot = 290kms-'). 
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Table 5. Model Results: Elemental Abundances. 
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Figure 10. Present epoch radial distributions of the logarithmic surface density of the atomic gas for 4 different mass distributions following the Vrot's values 
from each panel. In each one of them 10 models are represented, for the 10 different efficiencies, following labels in panel a). 



The consequence of a shorter collapse time-scale for the more 
massive spirals is clearly seen: the maximum is located at radii fur- 
ther away from the centre due to the exhaustion of the diffuse gas in 
the inner disc which moves the star formation outside. The smaller 
the galaxy mass, the closer to the centre is the maximum of the 
distribution, which resembles an exponential, except for the inner 
region. In fact, a shift in the maximum appears in each panel. In 
some cases, however, the low values of the surface gas density are 
due to the fact that the gas did not have enough time to fall com- 
pletely onto the equatorial disc. This effect is very clear for the 
model with A = 0.03 which shows a very steep distribution with 
gas densities lower than 15 M^/pc^. In this case the gas shows a 
radial distribution with a maximum at the centre. The low density 
out of this central region is not due to the creation of stars but to the 
long collapse time-scale which prevents the infall of the sufficient 
gas to be observed. 

Besides the variations due to the differences in total mass, 
which correspond to different collapse time-scales, and because 
we have selected different efficiencies to form stars and molecu- 
lar clouds, a same total mass may have produced discs in different 
evolutionary states. 

Thus, a same A = 0.15 may result in a disc of 7-8 kpc and 
atomic gas densities arotmd 5 MqIpc^, or a disc of only 5 kpc with 



a maximum density of 40 M^/pc^ in the region of 2 kpc. In the 
same way, a galaxy with a large value of the total mass, as the 
one with A = 2.50, may show a high gas mass density and a little 
disc for a high value of N or, on the contrary, be very evolved and 
therefore show no gas and a large stellar disc if the value of N is 
low. The first object (N > 5 ) could correspond to the case of low 
surface brightness galaxies, while the last ones (A' < 4) could be 
identified as the t5?pical high surface brightness spiral galaxies. 



3.2.2 Molecular gas radial distributions 

An important success of the multiphase models has been the abil- 
ity to reproduce the radial distributions for the atomic gas and the 
molecular gas separately, which is possible due to the assumed star 
formation prescription in two steps, thus allowing the formation of 
molecular clouds prior to the appearance of stars. 

Another important consequence of this SFR law is that it takes 
into account feedback mechanisms, even negative. If molecular 
clouds form before stars, this implies a delay in the time of star for- 
mation. The molecular gas shows an evolution similar to that of the 
diffuse gas, but with a certain delay. This delay allows the mainte- 
nance of a radial distribution with an exponential shape, as usually 
observed, for a longer time, although in some evolved galaxies H2 
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is also consumed in the most central regions, thus reproducing the 
so-called central hole of the molecular gas radial distribution, ob- 
serve d in some galaxies iNishivama & Nakaill200l| : lRegan et all 
l200lh . the MWG being one of them (see Fig6b in next Section). 

This is seen in Fig ll ll for efficiencies corresponding to A' > 5, 
(depending on the total mass), for which the model lines turn over 
at the inner disc, which corresponds to the regions located at the 
border between bulge and disc. Thus, the galaxies with the highest 
efficiencies (N < 5) show a maximum in their radial distribution of 
H2, which is always closer to the centre than that of the atomic gas 
distribution. Galaxies with the lowest efficiencies (A' > 5) show 
larger surface densities of molecular than atomic gas because the 
efficiency to form stars from molecular clouds is smaller than the 
efficiency to form these clouds. Therefore, the conversion of diffuse 
to molecular gas occurs more rapidly that the subsequent formation 
of stars. 



3.2.3 The radial stellar discs profiles 

The total mass converted into stars forms out the stellar disc in 
each galaxy. These stellar discs are represented in Fig |12l as the 
stellar surface density radial distributions. These profiles may be 
compared with brightness surface radial distributions, after the ap- 
propriate conversions. 



The values of efficiencies have a small influence in the result- 
ing stellar surface density distribution shape: the total mass of stars 
created is similar for all N < 6, although they are formed at dif- 
ferent rates, that is the resulting stellar populations have different 
mean ages. For the most evolved cases, most stars were created 
very rapidly, while for the less evolved ones, stars formed later as 
average, as it is shown in Fig.|7|for the characteristic radius region. 
Therefore, the radial distributions of surface brightness would re- 
sult very similar for galaxies of a given galactic total mass, but 
colors are expected to be different, redder for the galaxies with the 
highest efficiencies for the formation of stars and molecular clouds. 

A very interesting result is that the central value of the dis- 
tribution is practically the same, around 100 M^/pc^, for all rota- 
tion curves and efficiencies, in agreement with Freeman's law. Only 
galaxies with the smallest efficiencies or the less massive discs 
show central densities smaller than this value. We cannot compute 
a surface brightness only with these models, but assuming a ra- 
tio M/L = 1 for the stellar populations, this implies a surface lu- 
minosity density of ~ lOOLg/pc- and scale lengths in agreement 
with observed generic trends. Only in models for N> 7 the stellar 
discs show a different appearance: they look less massive, as cor- 
responding to discs in the process of formation. This implies that 
the surface brightness is lower than for the other types for a similar 
characteristic total mass. 
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Figure 12. Present epoch radial distributions of total mass surface density for 4 different mass distributions. Symbols meaning in panel a). 



In any case, we remind that all information related to pho- 
tometric quantities, and this also applies to the disc scale lengths, 
which may be obtained from the star formation histories and en- 
richment relations of models, need the application of evolutionary 
synthesis models, what is out of the scope of this work. 



3.2.4 The elemental abundances 

One of the most important results of this grid of models refers to 
the oxygen abundances, shown in Fig. ll3l A radial gradient appears 
for most of models. This is due to the different evolutionary rates 
along radius: the inner regions evolve more rapidly that the outer 
ones, thus steepening the radial gradient very soon for most model 
galaxies. Then, the radial gradient flattens for the more massive 
and/or most evolved (small N) galaxies due to the rapid evolution, 
even in the outer regions, which produces a large quantity of ele- 
ments, with the oxygen abundance reaching a saturation level. This 
level is found to be around 12 + log{0/H) ~ 9.0 - 9.1 dex. Ob- 
servations in the inner disc of our Galaxy support this statement 
( Smartt et al., 2001). This result was already found in our previous 
works. 

Moreover, the larger the mass of the galaxy, the faster the ef- 
fect: a galaxy with V,o, = 100 km.s"' has a flat radial gradient 
for efficiencies corresponding to N= 1 or 2, while a galaxy with 



Vror = 200 km.s"', shows a flat distribution for A' < 4. The less 
massive galaxies maintain a steeper radial distribution of oxygen 
for almost all efficiencies, with very similar values of the gradients. 

Nevertheless, for any galaxy mass, if A' > 7-8 the radial abun- 
dances distribution are flat. Thus, the less ev olved galaxies show no 
gradient, such it is observed in LSB galaxies Blok & McGaughL 
[lj)96J), the intermediate ones show steep gradients, and the most 
evolved galaxies have, once again, flat abundance radial distribu- 
tions. The largest values of radial gradients correspond to the in- 
termediate evolutionary type galaxies, with the limiting N varying 
according to the total mass of the galaxy. The more massive galax- 
ies only show a significant radial gradient if N< 8 while the less 
massive ones have a flat gradient only if N 5^ 9 with the rest having 
very pronounced radial gradients even for N = I. 

For the evolved galaxies, the characteristic efficiencies are 
high for all the disc, thus producing a high and early star forma- 
tion in all radial regions. In this case, the oxygen abundance reaches 
very soon a saturation level, flattening the radial gradient developed 
at early times of the evolution. The characteristic oxygen abun- 
dance, measured at R^, is higher for the more massive galaxies and 
lower for the less massive ones. However, this correlation is not 
apparent when the central abundance is used, due to the existence 
of the saturation level in the oxygen abundance, which produces a 
flattening of the radial gradient in the inner disc, even for galaxies 
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Figure 13. Present epoch radial distributions of Oxygen abundance, 12 + Iog{0/H), for 4 different mass distributions. Symbols meaning in panel a). 



with intermediate efficiencies. Actually, this saturation level would 
correspond to the integrated yield of oxygen for a single stellar pop- 
ulation. As oxygen is ejected by the most massive stars of a genera- 
tion of stars, it appears rapidly in the ISM. It is therefore impossible 
to surpass this level of abundance. 

In fact, the oxygen abundance radial distribution shows some- 
times a bad fit to observations in the central parts of the discs, 
which give frequently abundances larger than 9.10 dex. This ab- 
solute value of the oxygen abundance is not reached in any case by 
the models. All the computations performed within the multiphase 
approach reach a maximum 12 + log(0/H) ~ 9.10 dex which no 
model can exceed. We might think that uncertainties in the calcu- 
lations of the stellar yields elements, still very dependent on the 
assumed evolution of stars, are a possible reason for this saturation 
level appears in theoretical models. However, low and intermediate 
mass stars, for w hich yields show la rge variations among the dif- 
ferent groups (see lGavilan et alil2004i. and references therein), de- 
pending on the assimied hypothesis in the calculation of the stellar 
evolution in the latest stellar evolutionary phases, do not produce 
oxygen; these calculations only affect to N and C abundances. On 
the other hand, the most important uncertainty for the production 
of oxygen in massive stars refers to the strength of stellar winds 
responsible for stellar mass loss. This mass loss implies less pro- 
duction of oxygen and a larger ejection of carbon. In our models we 



have used Woosley & Weaver's yields which do no include stellar 
winds. Therefore, from this point of view, the modelled oxygen 
abundances must be considered upper limits. 

It should be recalled, however, that all oxygen data yield- 
ing values larger than 9.1 dex have been obtained from observa- 
tions of Hi: regions where the electronic temperature could not 
be measured, and hence the oxygen abundances have been de- 
rived through empirical calibrations which are very uncertain in 
the high abundance regime. Actually, the shape of the radial dis- 
tribution of oxygen changes depending on the calibration used 
(e.g. Kennicutt & Garnett . 1996). The suspicion that these abun- 
dances are overestima ted at least by 0.2 dex is very reasonable 
( |Pilvuginl l200ol l200lh . Abundance estimates in Hn regions his- 
torically considered as metal-rich, result to be almost solar once 



the ele ctronic temperatur e has b e en finally meas ured ( Diaz et al 
200Qt ICastellanos et all I200S iKennicutt. Bresolin. & Garnett . 



200alGarnett. Kennicutt. & Bresolirll2004 . 



3.2.5 The star formation rate 

Radial distributions of star formation rate surface density show an 
exponential shape in the outer disc, but a less clear one in the inner 
regions, where some models show a distribution flatter than that of 
the molecular gas, and even decreasing toward low values of the 
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Figure 15. The Solar Neighborhood evolution as results from the chosen 
JVIWG model for the region located at S = 8 kpc. a) The star formation 
history with data from Twaroa 1 1980) -filled dots- and iRocha-Pinto et al. 
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star formation rate in the centre. These SFR radial distributio ns are 
very similar to those estimated by 'Martin & KennicutJ i200lh from 
Ho- fluxes, as we show in Fig^|for a few galaxies. 

The star formation rate assumed in our models do not produce 
bursts of massive stars in the low mass galaxies whatever the effi- 
ciencies. Only the massive galaxies are able to keep a large quantity 
of gas in a small region, usually at the centre although sometimes at 
the inner disc regions. On the contrary, low mass galaxies collapse 
very slowly, and thus the star formation rate maintains a low level 
during the whole life of the galaxy. In fact, recent works suggest 
the same scenario for both low ma ss and lo w surface brightness 
galaxies ivan den Hoek et alH2000l : lLegrandL 12000: Braun, 200l|) 
in order to take into account the observed data. Our resulting abun- 
dances and gas fractions for lo w mass galaxies seem to b e in rough 
agreement with these findings iOavilan & Moll a.'2004b\ although 
our model results can not still be compared with photometric data, 
for which many more observational sets exist. 



and ll6l together with the available data. Fig. llSl shows the time evo- 
lution of the star formation history, panel a), and the age-metallicity 
relation, panel b), of the region located at R=8 kpc from the Galac- 
tic centre, as compared to data for the Solar Neighbourhood. The 
observational trends are reproduced adequately, although the inodel 
predicted maximum of the star formation rate appears slightly dis- 
placed toward earlier times with respect to observations. 

Fig. ll6l shows the present time radial distributions of diffuse 
and molecular gas, mass of stars, star formation rate and oxygen 
and nitrogen abundances for the galactic disc. The diffuse gas ra- 
dial distribution, panel a), is well reproduced in shape, although 
the maximum observed gas surface density is somewhat displaced 
to outer radii compared with the ol dest data. It, h owever, fits well 
the most recent data obtained from lNakanishi & Sofue 1.2004,) for 
radii R> 5 kpc. Regarding the molecular gas density, panel b), the 
distribution is quasi-exponential from ^ 8 kpc and decreases at 
the inner disc regions. Taking into account that recent data give 
low densities at these inner regions, we consider that our model 
results may be adequate. The stellar mass distribution, panel c), is 
exponential in shape in agreement with the surface brightness dis- 
tribution, and radial distributions, panels e) and f), reproduce those 
shown by the most recent data. 

The only feature which is not well reproduced by the model 
is the SFR radial distribution shown in panel d), which decreases 
toward the inner disc in apparent discrepancy with observations. 
The modelled star formation rate distribution has a maximum in 
/? ~ 7 - 8 kpc, while the observed one increases exponentially 
toward the galactic centre, or levels off at 3-4 kpc. On the other 
hand, the maximum of the atomic gas density is observed around 
10-11 kpc, and the molecular gas density has its maximum ob- 
served at around 6 kpc. Therefore, it results difficult to explain, 
from an observational point of view, how the star formation rate 
remains so high at the inner disc (inside the central 3-4 kpc), where 
both gas phases are already consumed. In fact, the recent data from 
fwilliains & McKee 1 1997) show a decline for R< 5 kpc and a de- 
creasing star formation rate for the inner disc regions is also ob- 
served in a large number of galaxies, as we will see in the next 
section. Actually, the MWG radial distributio n of the present star 
formation rate is still a matter of discussion ^S trong et ali |2004|) 
since the detection of young sources, which are embedded in gas 
clouds, is difficult and might be affected by selec tion effects. Data 
from pulsars, supernovae or OB star formation iBronfman et ali 
l200d) seem to indicate that the radial distribution of the SFR has a 
maximum around 5 kpc, decreasing toward both smaller and larger 
radii. On the other hand Case & Bhattacharva 1 1996) obtained a 
much ffatter distribution from the observed cosmic rays and gamma 
radiation spectra. Since a decrease at the inner disc is not unreason- 
able from the comparison with other spiral galaxy data, we consider 
that the model results for the MWG are acceptable. 



4 DISCUSSION 

4.1 Calibration with the MWG 

The first application of any theoretical model is to check its validity 
for the MWG. A large set of observational data for the Solar Neigh- 
bourhood and the galactic disc exists and therefore the number of 
constraints is large compared to the number of free parameters of 
the computed models. The model for N = 4, and total mass distri- 
bution number 28, corresponding to A = 1.0 and maximum rotation 
velocity V„,„j- = 200 km s"', is the model more representative of the 
MWG. 

The results corresponding to this model are shown in Figs. 1151 



4.2 Comparison with individual galaxies 

There is only a small sample of galaxies for which large observa- 
tional data sets, including neutral and molecular gas distributions, 
exist. In what follows we show a comparison of our model results 
with the data corresponding to theses galaxies. 

The characteristics and corresponding input parameters of this 
galaxy sample are given in Table|6| For each galaxy. Column (1), 
the morphological type index is given in Column (2), while the 
classical Hubble type is given in Column (3). The adopted distance 
(taken from references following TableQ Column (2) is in Column 
(4); the maximum rotation velocity is given in Column (5). The 
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Figure 14. Present epoch radial distributions of tlie suiface density of the star formation rate for 4 different mass distributions Symbols meaning in panel a). 
Table 6. Galaxy Sample Characteristics and Model Input Parameters. 



Galaxy 


T 


Type 


D 


^sliot.max 


Mass Distr. 


Tc 




N 


Told 






Name 




Class 


(Mpc) 


(kms"') 


Number 


(Gyr) 


(kpc) 




(Gyr) 






NGC 300 


7 


Scd/Sd 


1.65 


85 


13 


13.31 


2.3 


7 


13.3 


0.07 


0.007 


NGC 598 


6 


Sc/Scd 


0.84 


110 


21 


8.13 


2.9 


6 


10.3 


0.05 


0.005 


NGC 628 


5 


Sc 


11.4 


220 


31 


3.43 


7.7 


5 


3.28 


0.25 


0.01 


NGC 4535 


5 




16.6 


210 


30 


3.59 


7.4 


5 


3.50 


0.28 


0.02 


NGC 6946 


6 


Scd 


7 


180 


25 


4.94 


6.2 


6 


4.26 


0.18 


0.02 



number of the radial distribution of total mass -corresponding to 
the column (1) of Table ^ chosen to represent each galaxy is in 
Columns (6). The characteristic collapse time and the characteristic 
radius corresponding to each distribution are given in columns (7) 
and (8). The number N chosen as the best one to reproduce the 
observations is in column (9). The last columns, (10) to (12) give, 
for comparison purposes, the collapse time-scale and efficiencies 
for molec ular cloud and star formation used in our previous models 
jMoUa et al., 1996. 1999) for these same galaxies. 

The radial distributions of the different quantities: atomic and 
molecular gas densities, star formation rate and oxygen abundance, 
for these galaxies are shown in Fig. 1171 together with the corre- 



sponding observational data, taken from the references given in Ta- 
bleQ 

In the first row of panels of Fig. ll7l we can clearly see that the 
radial distributions of neutral hydrogen are very well reproduced 
by the models. For each galaxy, the distribution shows a maximum 
along the disc, as predicted. Also, for galaxies with similar total 
mass, such as NGC 4535 and NGC 6946, this maximum is higher 
for later morphological type. For galaxies with the same value of 
T, but different total mass, such as NGC 6946 and NGC 598, this 
maximum does not change its absolute value, but it is located at 
a different galactocentric distance, closer to the galaxy centre for 
the less massive one, due to the longer collapse time-scale which 
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Figure 16. Present epoch radial distributions for the MWG simulated galaxy (distribution number 28, A = 1.00, T = 4): a) atomic gas density with data from 
[Garwood & Dickey 1 1 989 ) -GW-, Wouterlook et al. (1 990 ) -WOU-, Rana 1 1 99 1 ) -RA- and Nakanishi & Sofue 1 2004 ) -NAK; b) molecular gas surface den- 
sity with data from Grabelsky et al. 1 1987) -GR, Bronfman et al. 1 1988) -BR-, Wouterlook et al. 1 1990) -WOU- and Williams & McKee ( 1997) -WILL-; c) 
stellar surface density with data from Talbot 1 .198 0); the dashed line corresponds to .Vallen ari^ertell^& Schmidtobreic^ i^QO^); d) star formation rate surface 
densit y normalized to the present time solar value in logarithmic scale, data taken from lLacev &'Fa!Tn98^0GUstei^feN ^^^ 11^983); Williams & McKa" 



1^ 



19^), LF, GM and WILL, respectively; e) and f) oxygen and nitrogen abundance as 12 + log(X/H) with data from Afflerbach, Churchwell, & Werna 
) -AFF- , Esteban et al. 1 1999); Esteban, Peimbert, & Torres-Peimberl ( 1999b); Esteban. Peimbert. Torres-Peimbert. & Garci'a-Roias 1 1999c) -EST-, 

[Fich & Silkey ( 1 9911) -FS- . Fitzsimmons et at, |1992|1 -FZ-. iPeimbert |1979|) -PE-,|Shaver et al. 1 1983) -SH- , Vflchez & Esteban 1 1996), -VIES- and 
Smartt & RollestoiJ_^^9^^martt et al. ( 2001) -SM- as lab elled in panel e) . Large fille d symbols at R= 8 kpc in both panels represent the sola r abundance s 
from lGrevesse & Sauv 3IT998) - circles-, Holwegei 1 2001) -squares- and Allende Prieto. Lambert & Asplund ( 2001 ); AUende Prieto et al. ( 2001 , 2003), 
crosses- ; while th e large empty sy mbols are the inters tellar medium abundances given b y Meve r. Cardelli & Sofia 1 1997) ; Mever. lura & Cardelli i , 1998i) . 



-crosses- ; wtule tn e large empty sy mbols are tne inters tellar medium abunaances given b y 
-circles-. |p'eimberlllT999l) -sauares-. ISofia et al]<200lil . -triangles, and lMoos et"aiTi2002l) . 
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Figure 17. Present epoch radial distributions for atomic and molecular gas densities, in units of Mo/yr, (first and second rows), star formation rate in 
MoGyr'pc"-, (third row), and oxygen abundance 12 + log(0/H) (last row) for the sample used to check the grid of chemical evolution models (Tablel^. The 
observational data are taken from references given in TablelTl 



reflects on a slower evolution. We would like to point out that the 
agreement between model results and data is much improved when 
good quality data (usually the most recently published ones) are 
used. The selection of the best obtained distances improves extraor- 
dinarily this agreement, as well. The same applies to the rest of the 
panels in this figure. 

The radial distributions of molecular cloud surface density for 
each galaxy, except NGC 300 for which no data exist, are shown 
in the second row of Fig. 1171 The agreement between model results 
and observations is good for the outer discs which follow a quasi 
exponential distribution, but is not so good for the inner discs where 



the molecular hydrogen surface density is observed to continue in- 
creasing instead of turn over as predicted by the models. We might 
artificially decrease the efficiencies en in these zones, instead to 
maintain them constant as we do, and thus recover the observed ex- 
ponential shape. But in this case oxygen abundances will be smaller 
than observed in these same regions. On the other hand we should 
recall that the molecular masses are estimated from the CO inten- 
sity through a calibration factor which depends on metallicity in 
a way which would produce smaller molecular gas densiti es than 
usually assum ed at the inner galactic disc ( Verier 8l HodgtA ll995t 
IWilsor[[l993) . Taking into account that recent data yield low den- 
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sities at these inner regions, we consider that our model results 
may represent adequately the reality. In any case, differences be- 
tween models and data are larger than in the case of the diffuse gas, 
which is not unexpected, given the larger uncertainties involved in 
the derivation of molecular hydrogen masses. 

The radial distribution of the star formation rate for each 
galaxy is shown in the third row of Fie ll7l The agreement between 
model results and data is very good for the most massive galaxies 
for which both the maximum of the star formation rate and its loca- 
tion is well reproduced. For the less massive galaxies, the predicted 
central turnover of the distribution is not observed. 

In the last row of Fig. 1171 the oxygen abundance radial distri- 
bution for each galaxy, except NGC 4535 for which no data exist, 
is shown. The oxygen radial gradient is reproduced in all cases and 
also the observed trend of steeper radial distributions - larger ra- 
dial gradients- for the late type galaxies is well reproduced by the 
models. 

We would like to emphasize that the models shown in this 
section have not been computed specifically for each of the sample 
galaxies, as was done in our previous works. The models shown 
in Fig^|were selected among the 10 available ones for their cor- 
responding rotation velocity following Tabl^of the grid. Thus the 
good agreement found between model results and data demonstrate 
that our bi-parametric models are able to adequately reproduce real 
galaxies. 

4.3 Star formation rate vs gas surface density 

The computed star formation rate reproduces the relation obtained 
byjKennicutt ( 1989), when it is represented vs the total gas surface 
density as can be seen in Fig. 1181 In the two upper panels of the fig- 
ure (A^ ^ 6), we have over-plotted the data from Young et al..(.1996.) 
and lKennicutil<I99St) which are seen to fall in the locus defined by 
the models. Models with N > 7 seem to be out of the region where 
this data lie. However, due to their extremely low efficiencies, we 
may assume that these models do not simulate normal bright spi- 
ral galaxies, but other kind of objects with a low stellar content 
and high gas fractions more simi lar to dwarf and LSB galaxies. We 
have, therefore, taken data from lLeerandl fcOOO) on these kinds of 
objects and computed the densities for both quantities assuming an 
optical radius of 5 kpc for all of them. Obviously, a change in this 
radial dimension would vary the final values of our estimates, but 
our hypothesis is probably valid within a factor of 2 (radii less than 
10 kpc). Under this assumption, we see that the points, shown in the 
third panel of Fig^| fall in the upper locus defined by the models. 
If the radius were smaller than assumed, the densities would be 
even higher, and the points would move in the figure following the 
direction given by the arrow. A second factor not included in these 
estimates, is the molecular gas content. There are some works sug- 
gesting that the molecular gas amounts to less than 10% in this kind 
of galaxies while other seems to indicate that the molecular mass 
may be as large as in the brightest massive spirals. In any case, the 
inclusion of this factor would move the points to the right. In both 
cases the data points would populate the region of the diagram oc- 
cupied by the models. We therefore conclude that our models are 
able to reproduce the observed trend of SFR vs total gas surface 
density for different types of galaxies. 

4.4 Radial gradients of abundances 

The model computed abundance radial gradients for oxygen have 
been obtained by fitting a least squares straigh line to the oxygen 
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Figure 18. The relation of the surface density of the star formation rate 
with the total gas density for 4 different mass distributions: a) models with 
N < 3; b) models with 4 ^ W < 6; and c) models with N > 1; Data are 
from fyoung et al. 1 1996) -open squares- and Kennicutt 1 1998)-filled dots-, 
in panels a) and b), while those from panel c) are given bv .Legrand 1.2000.1 . 



abundances for radii 0.5Rckpc < R < 3Rc kpc. The chosen ra- 
dial range tries to eliminate the central region (where the oxygen 
abundance distribution flattens) and the outer region where there 
are no data. Thus the calculated gradients may be compared more 
precisely with the corresponding observed ones. 

These radial gradients, measured as dex/kpc, are represented 
as a function of rotation velocity in Fig. 1191 The relation between 
both quantities shows that the radial distributions are steeper for 
the models with lower rotation velocities if N< 7. The radial gra- 
dients for a given rotation velocity are larger in absolute value for 
increasing N, and tend to zero value for the most massive galaxies 
and low N. The models with N > 7, however, tend to deviate from 
this function and show lower absolute values. That is, for a simi- 
lar rotation curve, their radial distributions of abundances show a 
flattening as compared to those corresponding to the models with 
N ^7 . The rotation velocity at which the models deviate from the 
common locus depends on A': for N = 7 it occurs for Vopi ~ 100 
while for A' = 10 almost all rotation velocity curves produce flat 
radial gradients of abundances. 

Data from Zarit.skv et all <I994 : iGarnett et all ll997i) : 
Ivan Zee. Salzer & Havnes, (il99a) are over-plotted on this Fig. 1 191 
We see that the modelled trend reproduce the observations, al- 
though a more profound analysis about this (and other) correlation 
will be performed in the future. 



4.5 Efficiencies 

In order to use this grid for a given galaxy for which observational 
data are known, we must first select the radial distribution of mass. 
This means that we must know the total mass, the maximum rota- 
tion velocity or, if none of them is available, the luminosity or the 
magnitude in the I band. According to this value we may choose the 
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Table 7. Galaxy Sample Data References. 



Galaxy 



Hi 



H2 



SFR 



[0/H] 



NGC 300 
NGC 598 
NGC 628 
NGC 4535 



TU 
TP 
MK 
TP 



NGC 6946 WE 



ROG, PUC 
NEW, CS 

WEV 
CAY, WAR 
CAR 



HEY 
AL, NIS, RY, SS 
KY, NIS, RY 
RY, TY WAL 



DEH 
HEY, KEN, RK 
LR, MK, RY 

KEN, MK 
KEN, MK, RY 



PAG, CHR, DEH 

KA, VIL 
BR, MCC. ZEE 

BR, MCC 



AL: lAdi er & Liszl Il989t) : BK: IBelle^^^mJ 
phristensen e t alj 1997t): CA Y: ICayatte et aU f 



iBellev & Ro\< Il992l) : CAR: IC:arig nan et aJ 



1991 



1§2(J); CHR: 
ICorbelli & Salucci ( 20 0^; DEG: 
HEY: Hever et alj 



ai]U28 
lA: Ikw 



phi 

[oeeioia-Eastwood et al.' '1984'); DEH: Deharvena. Caolan. Leaue ux. et 

|M04); KY: Kennev & Youna il989); KEN: Kennicutt ( 198^; KA: |Kwitter & Allei ( 1981) ; 
LR: iLelievre & Rov (2000); MAR: Martin & Kennicutt (2001); MCC: McCall et alJ Il985h ; 
NN: iNishivama & Nakai (2001); NEW: Newton (1980); PAG: Paael et al. (1979|); PUC: 
l^he. Carianan. & Bosma (1990); ROG: Roastadetal. (1979); RK: Rumstav & Kaufman (1983); 
RY: .Rownd & Youn a 1 1999); SS: Saae & Solomon ( 1989); TY: .Tacconi & Youna (.1989.); TU:,Tull\j 
ll988t|: TP:ITu11v & Pierce ( 2000); ZEE: van Zee. Salzer & Havnesli 19981): VIL:lvi1chez et alJ<1988i-| : 
WAL: IWalsh et alj< 20021) ; WAR: lWarmel3tl988i) ; WEV: iWevers. van der Kmit & AlletJ h98a) 
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Figure 19. Radial gradients in dex kpc"' as obtained from the modeled oxy- 
gen abundance radial distributions, computed with a least square straigh line 
in the radial range 0.5Rckpc < R < 3Rckpc, as a function of the rotation 
velocity Vmax in kpc.s"'. Models are represented by the same line cod- 
ing thari Fig. 10 to 14. Data are fr om Zaritsk v et alj [l994l) ; lGamett et al] 
ll997t) ; lvan Zee. Salzer & Havnej tl998) . shown by squares, dots and tri- 
angles, respectively. 



number of the distribution from TableQ Then, 10 different models, 
corresponding to the 10 different efficiencies, are available. 

The standard procedure in chemical evolution would be to 
see which of them is able to reproduce with success known data 
such as elemental abundances or gas densities. If more than one 
observational constraint is available, some kind of minimum error 
or maximum probability technique may be used for the purpose 
of choosing the best model. Once this selection is performed we 
may use the time evolution given by the chosen model to predict 
the star formation history, the age-metallicity relation, the stellar 
populations or any other quantity relative to the modelled galaxy. 



Obviously, the largest the observational number of constraints, the 
smallest the uncertainty in the selection of the best model. The old 
and well known uniqueness problem of the chemical evolution sug- 
gest that more than one model may reproduce the data. Actually, 
we have shown in MoUa & Hardv ( 2002 ) that this problem reduces 
greatly when more than two observational constraints can be used. 
Only 4% of the 500 models computed with with different input pa- 
rameters could reproduce the observations relative to HI, oxygen 
abundance and SFR at the same time. Therefore we are confident 
that models may provide the evolutionary history of a given galaxy 
within a reasonable accuracy. 

Sometimes, however, these observational constraints are not 
available (or not with sufficient precision) to select the adequate 
model. In that case, some other method to choose the possible evo- 
lutionary track is necessary. Our best models from previous works 
give us evidence that the efiiciencies to form molecular clouds 
from atomic gas and stars from molecular gas, seem to depend on 
the galaxy morphological type. If this were the case, the selection 
of the best model might be easie r. In f a ct, this dependence wa s 
already found in Ferrini & Galll Jl988h : lOalli & Feriinil Jl989h . 
where these authors quantified the efficiencies to form molecular 
clouds and the frequency of cloud-cloud collisions, finding that a 
variation of 10 in the parameters H and /j is needed when the Hub- 
ble type changes from one stage to the next (Sa to Sb, etc.). 

In order to check if a relation between efficiencies and mor- 
phological type exists, we have plotted in Fig. 1201 panel a) the log- 
arithmic efficiency e^, computed following the equation (19) from 
section 2.3, as a function of T for the data of Young et al. ( 199^. 
There + symbols correspond to the raw values while solid dots rep- 
resent the averaged values obtained for 10 bins, one for each T. A 
clear decreasing correlation appears for T < 9. Values for T > 9 
fall above the trend but the number of points is small for what we 
assume that the final increase is probably spurious. Based on our 
previous works, we expect a relation of the form: oc exp"^^^^. 
Therefore, in panel b) we have plotted the logarithmic efficiency 
as a function of where a linear correlation is apparent. A least- 
squares fit gives, for T <9 : 



< Ine,, >= 



- 1.60-0.0577 



1.19 



Therefore: 



(23) 



(24) 
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Figure 20. a) The efficiency Me^ as a function of the morphological type 
T represented by + for all galaxies. The solid dots are the averaged val- 
ues computed with 10 bins. The solid line is a least-square second degree 
polynomial, b) The efficiency /ne^ as a function of T^. The solid line is the 
least-squares straight line. 

On the other hand, taking into account the relationship ob- 
tained for the ratio between both efficiencies, e,, and e^, this latter 
one may also be expressed by a similar function : e/j oc exp"^''^ 
where B would be 8, following our hypothesis from Section 2. 

The result that efficiencies are related to galaxy morphological 
type may be useful when a model for a given galaxy must be chosen 
and the observational data are scarce. In that case, once the total 
mass distribution has been selected from the rotation curve or the 
luminosity for the galaxy, an initial choice can be made under the 
assumption that the galaxy evolution is represented by the model 
of efficiencies N with N = T. We warn, however, that given the 
large dispersion of the data around the averaged efficiency values 
in Fig. 1201 other values of N are also possible. 



This finding would imply that the features of spiral galax- 
ies depend at least on two parameters: total mass and morpho- 
logical type, as some st udies about the characteri stics of galaxies. 



as the seminal work by 


Roberts & HavnesI Jl994). show (see also 


Vila-Costas & Edmunds 


,ll992llZaritskv et allll994lGamett et all 



1997, among those more related with the chemical evolution of 
galaxies). 



It is also well known that the radial gradient, measured as 
dex kpc ', observed in spiral galaxies depends on morphological 
type: late-type galaxies show steeper radial distributions of oxygen 
abundances than earlier ones, which showjnjomecases almost no 
gradient iOev & Kermicut t. 1993; D util & Rovlfl999n . At the same 
time, a correlation seems to exist between oxygen abundance and 
total mass surface density which, in turn, is related to morpholog- 
ical type. This correlation is stronger when the total mass of the 
galaxy, including the bulge, instead of the mass of the exponential 
disc alone, is used ( Rvder, 1995). It is evident from these works that 
a linear sequence with the mass can not be obtained, and therefore 
our models, being a bi-parametric grid, would be adequate to find 
the best model for each galaxy, although the existing relation be- 



tween the morphological type of galaxies and their total mass does 
difficult to discriminate which of these features is the origin of the 
different evolution of galaxies. 

From the chemical evolution point of view, a recent work 
(MoUa & Marquez, in preparation), where the multiphase evolu- 
tion model has been applied to a large sample (67) of spiral galax- 
ies, seems to demonstrate that the total mass, through its influence 
on the collapse time scale, has a larger effect on the evolution of a 
galaxy and its final radial gradient of abundances than the selected 
values of the named efficiencies. Thus, it seems as if galaxies would 
evolve mostly due to their total mass with some dispersion around 
the mean trend depending on hydrodynamical and environmental 
characteristics, which are taken into account in some way by the 
efficiencies values. The morphological type would then be the con- 
sequence of the different conditions of the intergalactic medium out 
of which the galaxy formed. 



5 SUMMARY AND CONCLUSIONS 

We have calculated the chemical evolution of a wide set of theoret- 
ical galaxies characterised by their total mass, through the collapse 
time scale, and their efficiencies to form molecular clouds and mas- 
sive stars from cloud-to-cloud collisions. 

With the selection of parameters and inputs described above, 
we have ran a total of 440 models, with 44 different rotation curves 
-implying 44 values of total mass, characteristic collapse time- 
scale and disc radius- and 10 sets of efficiencies for each one of 
them, implying 10 evolutionary rates for the star formation and gas 
consumption in the disc. 

For each model we have obtained the time evolution of the 
halo and the disc, and therefore the corresponding radial distribu- 
tions for the relevant quantities (masses, abundances, star formation 
rate, etc.). The star formation history for each radial region (halo 
and disc, separately) and within each one, the mass in each phase 
of matter: diffuse gas, molecular gas, low-mass stars, massive stars 
and remnants, has been followed. Besides that, we have obtained 
the abundances of 15 elements: H, D, He3, He4, C12, C13, N14, 
016, Ne, Mg, Si, Ca, S, Fe and neutron-rich nuclei, in all radial 
regions for both halo and disc. 

The results of our work can be summarised as follows: 

(i) The atomic gas shows a maximum in its radial distribution 
for all galaxies. This maximum is nearer to the galaxy centre in 
the low mass or less evolved galaxies than in the more evolved or 
massive galaxies, for which the maximum is along the disc and 
moving toward the outer zones. This behavior produces, in some 
cases, a central iioJe in the distribution of the diffuse gas. The dif- 
fuse gas radial distribution results to be a very strong constraint for 
selecting the best model out of the 10 computed ones with differ- 
ent efficiencies, corresponding to the total mass of a given spiral 
galaxy. 

(ii) The oxygen abundance reaches a maximum level, as a con- 
sequence of a saturation effect which occurs earlier for the mas- 
sive and more evolved galaxies. The less evolved galaxies do not 
reach this saturation level, except in the central region, and there- 
fore show a steep radial gradient in their oxygen abundance. The 
less massive and less evolved galaxies have not yet developed a ra- 
dial gradient and show flat radial distributions. This simulates an 
on-ofF effect: for A? = 7 a radial gradient appears if A > 0.15 while 
at W = 8 it only appears for A ~ 1.50. This behavior is in agreement 
with observations and solves the apparent inconsistency shown by 
trends showing steep gradients for late type galaxies and flatter ones 
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for the earliest ones, while, at the same time, most irregulars show 
no gradient at all and very uniform abundances. 

(iii) The model calculated star formation rates reproduce the ob- 
served trend between the surface density of this quantity and the 
total gas density including, not only the massive normal galaxies 
data, but also the low surface brightness ones. 

Actually, to our knowledge, there are not other chemical evo- 
lution models which compare predicted with observed radial distri- 
butions of diffuse and molecular gas, star formation rate and abim- 
dances for disc galaxies other than MWG. In this sense our models 
should be considered as an improvement over the standard ones. 

A study of the possible correlations among galaxy properties 
and different features obtained from the complete set of results is 
now possible, since a statistically significant number of theoretical 
models is available. This requires a deeper analysis which is out of 
the scope of this paper and will be done in the next future. 
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Table 5. Model Results: Elemental Abundances. 
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